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ABSTRACT 

Using the 100-m Green Bank Telescope, we have conducted a cm-wavelength 
I search for CO J=l-0 line emission towards the high-redshift, far- infrared luminous 

r~| . object, HDF850.1 over the redshift interval 3.3 < z < 5.4. Despite the wealth of existing 

^ i' multi-wavelength observations, and the recent identification of a galaxy counterpart 

in deep K' band (2.2 /im) imaging, an unambiguous spectroscopic redshift has not yet 
been obtained for this object. A far-infrared-to-radio wavelength photometric redshift 
technique however, predicts a ~90% probability that the redshift is in the range, 
I 3.3 ^ z <, 5.4 (equivalent to an observed redshifted CO J=l-0 emission line frequency, 

26.5 > i^obs ^ 18.0 GHz), making HDF850.1 a potential occupent of the 'high-redshift 
tail' of submm selected galaxies. We have also conducted a search for CO J=2-l line 
K> , emission over the narrower redshift range, 3.9 <, z <, 4.3. Although wc do not detect 

^ • any CO line emission in this object, our limits to the CO line luminosity are in broad 

I agreement with the median value measured in the current sample of high-redshift, 

submm selected objects detected in high- J CO line emission, but not sufhcient to fully 
test the validity of the photometric redshift technique. 

Key virords: galaxies: starburst - galaxies: individual: HDF850.1 - radio lines: galax- 
ies " cosmology: observations 
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1 INTRODUCTION 

In recent years, blank-field extragalactic surveys at submil- 
limetre/millimetre (hereafter submm) wavelengths have re- 
vealed a population of dusty galaxies undergoing vigorous 
star formation in the young Universe (e.g. Smail, Ivison & 
Blain 1997; Hughes et al. 1998; Barger et al. 1998; Bertoldi 
et al. 2000; Scott et al. 2002; Borys et al. 2003; Greve et al. 
2004; Laurent et al. 2005). Independent methods of redshift 
determination for these objects imply that the majority lie 
at, z ;> 2 (Aretxaga et al. 2003, 2005; Chapman et al. 2003a, 
2005). Thus their inferred star formation rates are in the 
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range, 100 to 1000 M© yr~^, and arguably require large 
reservoirs of molecular gas for fuelling such high levels of 
sustained activity. Measuring the total molecular gas mass 
contained within the interstellar medium (ISM) of a high- 
redshift submm galaxy (hereafter SMG) is most effectively 
accomplished through observations of redshifted molecular 
CO line emission. 

Over the past 15 years, observing CO line emission in 
high-redshift objects has become a powerful means of con- 
straining the physical conditions of the gas within their 
molecular ISM (for an excellent review see Solomon & Van- 
den Bout 2005). The luminosity in the CO J=l-0 line is the 
optimal estimator of the total molecular gas mass yet, for 
practical reasons, in gas-rich objects at high-redshift (z ^ 2), 
searches are normally first conducted for emission from high- 
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J ( Juppcr ^ 2) CO line transitions, and subsequent searches 
for CO J— 1-0 emission are carried out if the mm- wavelength 
searches for high- J CO lines are successful. This approach 
can bias the sample of objects detected in CO line emission 
to those with hotter, and denser gas. 

To date, high- J CO line emission has been detected in 
14 SMGs (Frayer et al. 1998, 1999; Neri et al. 2003; Genzel et 
al. 2003; Downes & Solomon 2003; Sheth et al. 2004; Kneib 
et al. 2005 Greve et al. 2005; Tacconi et al. 2006), while 
the CO J=l-0 line has been detected in only one of these 
(Hainline et al. 2006). In addition to the faintness of molecu- 
lar emission lines in high-redshift SMGs, searches have been 
hindered by the limited spectral bandwidth of current mm- 
wavelength facilities, generally covering ~1700 km s"^ at 
3 mm, while the typical SMC CO linewidth is ~800 km s~^ 
FWHM. This narrow bandwidth can also be restrictive as 
galactic outflows in many high-redshift SMGs may lead to 
velocity offsets between the redshifts derived from Lya and 
CO emission lines. In some SMGs this difference may be 
greater than 600 km s^^ (e.g. Greve et al. 2005), due pos- 
sibly to galactic outflows, or scattering of Lya photons by 
dust. As CO emission- line frequencies are not expected to 
be biased with respect to the systematic redshift, broadband 
spectroscopic searches for CO line-emission should be a pow- 
erful means of obtaining redshifts for the SMG population. 

Educated searches for mm-to-cm wavelength molecular 
CO line emission in luminous, dusty galaxies without red- 
shifts, will become feasible in the near future as wideband 
spectrometers are available on large mm-to-cm wavelength 
telescopes, for example the 100-m Green Bank Telescope 
(GBlQ; Jewell & Prestage 2004), or the 50-m Large Mil- 
limeter Telescope (LM1[f|). In order to obtain redshift es- 
timates, and to guide the frequency tunings of these spec- 
troscopic searches for molecular line emission from SMGs, 
some groups have developed photometric redshift techniques 
which exploit the far-infrared-to-radio wavelength correla- 
tion in star-forming galaxies (Helou et al. 1985), or adopt 
template far-infrared spectral energy distributions (SEDs) 
based on nearby galaxies (Carilli & Yun 1999, 2000; Dunne, 
Clements & Bales 2000a; Rengarajan & Takeuchi 2001; 
Hughes et al. 2002; Aretxaga et al. 2003; Wiklind 2003; Hunt 
& Maiolino 2005). This technique has the potential to pro- 
vide redshift estimates for large samples of SMGs with in- 
dividual accuracies, Sz ~ ±0.3, when photometric flux mea- 
surements of three or more far-infrared-to-radio wavelengths 
are available (Aretxaga et al. 2005). 

The GET is the only operational mm-to-cm wavelength 
telescope in the world with instruments that have both suf- 
ficient spectral line sensitivity and receiver bandwidth to 
conduct guided searches for CO line emission in SMGs at 
redshifts z z 0.9. This lower redshift limit is set by the 
current GET frequency limit of 60 GHz and the CO J=l- 
line rest frequency of 115.2712 GHz. Given this restric- 
tion, an excellent candidate for conducting a blind search 
for CO line emission is HDF850.1 (Hughes et al. 1998), one 



^ The Green Bank Telescope is a facility of the National Radio 
Astronomy Observatory, operated by Associated Universities, Inc. 
under a Cooperative Agreement with the National Science Foun- 
dation. 

^ http://www.lmtgtm.org 
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Figure 1. The most recent submm-to-radio-wavelength photo- 
metric redshift estimate for HDF850.1, plotted as a function of 
redshifted (115.2712 GHz) CO J=l-0 line frequency (solid line). 
The dashed line shows an example of the opacity across the K- 
band due to the atmospheric water vapour line at ~22 GHz. 



of the most well studied SMGs, and the brightest 850 /im 
source in the confusion limited JCMT/SCUEA survey of 
the northern Hubble Deep Field. Due partly to the extreme 
optical faintness (K ~ 23.5, I — K > 5.2) of the gravita- 
tionally lensed galaxy counterpart to HDF850.1 (Dunlop et 
al. 2004), the redshift of this object has proven elusive. A 
wealth of deep rest-frame far-infrared-to-radio wavelength 
observations of HDF850.1 provide the basis for a photomet- 
ric redshift z=4.1±0.5 (Yun & CariUi 2002; Aretxaga et al. 
2003). In principle, HDF850.1 presents an ideal target for 
the GET with which to test the accuracy of our photomet- 
ric redshift technique, and thus has motivated a GET search 
for CO J=l-0 and CO J=2-l line emission over the redshift 
interval, 3.3 ^ z ^ 5.4. 

We present the results of a GET search for CO J=l- 
line emission in HDF850.1 over the redshift interval, 
3.3 ^ z ^ 5.4, and a search for CO J=2-l line emission over 
the narrower redshift interval, 3.9 ^ z <, 4.3. Throughout 
this work, we adopt the following A-dominated cosmological 
parameters: Ho = 70 km s~^ Mpc~\ = 0.7, Q,m = 0.3 
(Spergel et al. 2003, 2006). 



2 OBSERVATIONS AND DATA REDUCTION 

The far-infrared-to-radio wavelength photometric redshift 
estimate of Aretxaga et al. (2003) implies an 86 to 90% 
probability that HDF850.1 has a redshift in the range, 
3.3 ^ 2 < 5.4. Over this redshift interval, the 115.2712 GHz 
CO J=l-0 line is redshifted into the 18.0 to 26.5 GHz fre- 
quency window of the K-band receiver on the GET. Mo- 
tivated by this prediction, we have obtained a complete K- 
band spectrum of HDF850.1 in order to search for redshifted 
CO J=l-0 line emission. 

Observations in nod mode were carried out with the 
GET K-band receiver during October 2004 and May 2005. 
The position center adopted for the on-source beam was 
that of the mm-wavelength counterpart detected with the 
Plateau de Eure Interferometer (PdEI) by Downes et al. 
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(1999). All of the K-band observations were conducted under 
reasonably dry conditions, with an average 22 GHz zenith 
opacity, r22GHz ~ 0.09. The nearby quasar 3C295 was used 
for pointing purposes, as well as baseline and flux calibration 
throughout the K-band observations. To correct for slowly 
varying, large-scale spectral baseline features, the observa- 
tions were made by alternately nodding two beams sepa- 
rated by 178.8", between the source and blank sky. The GBT 
spectrometer allows a maximum instantaneous frequency 
coverage of 800 MHz bandwidth in each of 4 independent 
quadrants. For the observations presented here, one pair of 
quadrants was used to measure a ~1.5 GHz wide spectra 
on the source, while the other pair measured blank sky in 
the off-beam. A total of 6 tunings (or sequences) were there- 
fore used to cover the entire K-band window. Each spectral 
channel was 0.39 MHz wide so that the velocity resolution 
varied from ~4.4 km s"""^ to ~6.5 km s""*^ across the band. 
A total of 28.2 hours of integration time was devoted to 
the HDF850.1 K-band spectrum. The time spent on each 
~1.5 GHz tuning sequence was varied to compensate for 
the increased opacity towards the center of the band, due 
to the 22 GHz atmospheric water vapour line. The goal was 
to obtain a spectrum with uniform noise across the K-band. 
Overheads such as pointing, focusing, acquisition of baseline 
calibration spectra, and follow-up of potential CO line de- 
tections, amounted to an additional factor of 2 increase in 
the observing time. 

In December 2005, a search was also conducted for 
CO J=2-l emission using the GBT Q-band receiver (40 to 
48 GHz), over the redshift interval, 3.91 <, z <, 4.25. As 
in the case of the K-band observations, the nod observing 
mode was adopted, while the spectrometer was set up in 
wide bandwidth, low resolution mode. The velocity resolu- 
tion varied from ~2.5 km s~^ to ~2.7 km s~^ across our Q- 
band spectrum. As there is a dearth of bright, compact cali- 
bration sources at these higher frequencies, the primary flux 
calibrator was 3C286, while the nearby quasar 11534-495 
(~17° separation) was used for pointing. 

Obtaining quality spectra with the GBT Q-band re- 
ceiver generally requires low wind speeds (^3 m s~^), and 
an extremely dry, stable atmosphere. Only on a single night 
in December 2005, was data obtained under just such con- 
ditions, with a median r44GHz ~ 0.1 and negligible wind 
speeds. Despite acceptable Q-band weather conditions dur- 
ing this, and possibly one other observing shift, only the left 
polarization Q-band spectra is included in our analysis, as 
the majority of the right polarization spectra suffer from a 
severe baseline ripple of unknown origin (see i]2.1[) . 

Both the K-band and Q-band spectra were reduced us- 
ing the (^fetirflfl data analysis package. For a series of consec- 
utive scans, a co-added spectrum is produced following the 
standard procedures described by Vanden Bout, Solomon & 
Maddalena (2004), which are only outlined here for com- 
pleteness. 

For simplicity, we will consider the spectrum in only 
one quadrant, and one polarization from each of the 
two beams in the discussion that follows. Let us refer to 
these spectra as B\[v) and B2{i^), where one of the two 
beams is always pointed 'ON' the source, while the 'OFF' 

^ http://gbtidl.sourceforge.net 



beam is observing blank sky separated from the source by 
178.8"in azimuth. A single scan is the time that one beam 
spends on-source before nodding to the off-source position, 
when the other beam is nodded onto the source. For our 
observations, a scan duration of 1 minute, during which one 
beam is on-source before nodding to the off-source position, 
was adopted. This scan duration is chosen so as to minimize 
the frequency dependent variation in the sky brightness 
temperature between successive scans, while also spending 
less time on overheads such as nodding. Assuming that 
B\{u) begins on-source at scan i, then a series of normalized 
spectra are produced by subtraction of the off-source scan 
from the on-source scan (note that from this point on, the 
u dependance will not be made explicit): 
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etc.... 

It is important that each normalized spectrum be created 
from the 'ON'-'OFF' scan pair that travel down the same 
signal path (ie. the signal from a single beam and polar- 
ization), so that any path-dependent sources of baseline 
structure may be properly subtracted. An average of these 
normalized spectra is then calculated, with each spectrum 
weighted by the inverse square of the system temperature, 
Tgyg. The units of this mean spectrum are then calibrated 
from [K] to [Jy] using the best-fit gain curves shown for 
the K-band observations in Figure [21 The gain curves were 
derived from observations of objects with known (and non- 
variable) radio flux densities (3C295 for the K-band spec- 
trum, and 3C286 for the Q-band spectrum). 

The left and right polarization spectra are analyzed sep- 
arately so that any potential detection of CO line emission 
would have to be confirmed in both polarizations. A large 
fraction of the data (~53%) were considered unusable due to 
various forms of spectral baseline irregularities and contam- 
ination which could not be removed reliably during the data 
reduction process (see ij2.ip . We attempted to use spectra of 
the bright pointing sources to correct the baseline shapes in 
these spectra (Vanden Bout, Solomon & Maddalena 2004), 
however we found that this did not improve our results. Af- 
ter removal of the poor quality data, the total on-source 
integration time devoted to the final K-band spectra is 12.3 
hours in the left polarization spectra, and 13.9 hours in the 
right polarization spectra (see Table [l]). 

2.1 Spectral Baselines 

A major obstacle faced by searches for faint, broad emission 
lines, is the presence of variations in the spectral baseline 
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Figure 2. The GBT K-band gain (in Jy/K) for the left and right polarizations, averaged over all observations of the calibrator, 3C295 
{solid lines). Also plotted as dashed lines are the bestfit linear regression curves to the data. These fits are then used to calibrate the 
combined HDF850.1 spectra. 



Table 1. Summary of K-band observations of HDF850.1. The 
total integration time, tint, is based on the data which contributes 
to the final, left (L) and right (R) polarizaion spectra. 



scan: 63 pol: L bank: 
scan: 250 pol: L bank: 



Seq.# 


1^1 - 1^2 


GHz] 


iint (L/R) [hrs] 


a (L/R) [mjy] 


Kl 


22.10 - 


23.62 


2.7 / 3.1 


0.58 / 0.69 


K2 


23.54 - 


25.06 


2.2 / 1.8 


0.51 / 0.61 


K3 


20.88 - 


22.40 


2.1 / 3.1 


0.50 / 0.52 


K4 


19.44 - 


20.96 


1.8 / 2.6 


0.60 / 0.46 


K5 


24.98 - 


26.50 


1.5 / 1.2 


0.98 / 0.85 


K6 


18.00 - 


19.52 


1.9 / 2.1 


0.77 / 0.72 



shape, with a characteristic scale similar to the expected CO 
linewidths. These features are not uncommon in mm/cm 
wavelength spectra obtained with single-dish instruments, 
and can easily be mistaken for detections of weak, broad 
emission lines. There are a number of instrumental, as well 
as atmospheric effects which may produce spectral baseline 
artefacts. Here, we summarize a few such artefacts which 
have been identified in our GBT (K and Q-band) spectra of 
HDF850.1: 

• Weather variations over short timescales (on order of 
the length of a single scan of 1 minute duration) may result 
in inaccurate subtraction of the atmospheric contribution 
to the system temperature across the band. This is particu- 
larly problematic at Q-band frequencies, where clouds pass- 
ing overhead may lead to rapid changes in the temperature 
of the atmosphere within a single beam. 

• There may be interference from sources along the path 
of the analog signal, originating at the receiver on the tele- 
scope and travelling to the spectrometer backend in the 
GBT control room. An example of this is the 50 MHz rip- 
ple which appears in the left polarization K-band spectra 
taken during May, 2005 (Figure [SI, and is due to tempera- 
ture variations in the equipment room where the spectrom- 
eter is located. These variations cause standing waves in the 




100 200 300 400 

Frequency scale [MHz] 

Figure 3. a) An example of the ~50 MHz ripple which infected 
much of the left polarization, K-band spectra during the May, 
2005 observing run. For comparison, we plot a spectrum of the 
same frequency tuning and polarization obtained on a different 
date, b) The power spectra of the two spectra, exemplifying the 
prominant spike on scales of ~50 MHz in one of the two. This 
excess in power is likely due to temperature variations in the 
spectrometer room, which can result in the creation of standing 
waves in some of the signal connectors. 



connectors which manifest themselves as ripples in the spec- 
tral baseline. 

• Resonances in the receiver feeds may cause a loss of 
power at certain frequencies, leading to emission or absorp- 
tion lines, sometimes referred to as 'suck-outs'. When calcu- 
lating (ON — OFF)/OFF from the on-source and off-source 
spectra, there will be a feature which is either in emission 
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or absorption, depending on whether the power loss is in 
the on, or the off-beam. These 'suck-outs' are apparent in 
the K-band receiver temperature curve shown in Figure|4l[a) 
(e.g. at frequencies of ~22.6 GHz and ~25.7 GHz), and may 
also be present in the Q-band receiver system (however at 
present, no high- resolution receiver temperature data are 
available for this receiver). 

• A high-frequency ripple severely affected the right po- 
larization Q-band spectra taken in December, 2005 (but may 
also be present in the left polarization at a lower amplitude), 
and is of an unknown origin. 

Data which was infected by any of these artefacts was not 
included when creating the final spectra. Thus ~53% of the 
original data were discarded. 



3 RESULTS AND DISCUSSION 

We do not find any evidence for CO J=l-0 line emission 
in our K-band spectra. Figure |4ja) shows the final K-band 
spectra of HDF850.1 over the observed frequency range, 
Vohs = 18.0 to 26.5 GHz. Due to the presence of residual 
baseline features, ~53% of the raw data are not included 
in the final, co-added spectra shown here. The most notice- 
able contaminant of our spectra are the receiver resonance 
features, appearing at various frequencies along the K-band 
receiver temperature spectra (top curve, Figure|4ja)). These 
resonances originate within the feed horns and their ampli- 
tudes depend strongly on the weather conditions. Although 
these features prevent us from placing any constraints on the 
presence of CO J=l-0 line emission over certain redshift in- 
tervals, we are still able to obtain CO J=l-0 line luminosity 
limits over much of the K-band window. 

3.1 CO line luminosity limits 

We calculate 3-(t upper limits to the CO J=l-0 and CO J=2- 
1 line luminosities across the available K-band and Q-band 
spectra, respectively. We assume a range of CO line widths, 
AViinc ~ 460, 780 and 1100 km s~^, which represent the 
first, median and third quartiles of the linewidths in the 
first 12 SMGs detected in high- J CO line emission (see 
Greve et al. 2005) . The 3-cr upper limit to the CO line inte- 
grated intensity is given by (e.g. Isaak, Chandler & Carilli 

2004), 3- a/ AVlino/AVchannol ' ""channel ' AVchannol (in UUitS of 

Jy km s~^), where the velocity width of a channel, A14hannoi, 
and the r.m.s. per channel, (Tchannei, both vary across the 
combined spectra. The limits to the integrated line intensity 
are converted to 3-a upper limits to the CO line luminosity, 
L'qq, following the expression given by Solomon, Downes & 
Radford (1992). These limits are shown in Figure |4jb) for 
the CO J=l-0 line in the K-band spectra, and in Figure[5jb) 
for the CO J=2-l line in the Q-band spectrum. For compar- 
ison, we also plot the median CO line luminosity measured 
in the first 12 SMGs detected in high- J ( Juppor > 2) CO line 
emission (see Greve et al. 2005). 

The fact that neither CO J=l-0 , or CO J=2-l line 
emission is detected in HDF850.1, can be explained by three 
possible scenarios; i) the total molecular gas mass in this ob- 
ject is low, resulting in CO line emission that is weaker than 
our line luminosity limits, it) the frequency of the emission 



line coincides with that of a receiver resonance feature, or 
lii) the redshift of this object is such that the emission fre- 
quency of the CO J=l-0 line is outside the range accessible 
to the K-band receiver, a possibility which has a 10-14% 
probability according to the photometric redshift estimate 
of Aretxaga et al. (2003). 

Under the assumption that is the correct scenario, 
the limits to the CO J=l-0 line luminosity can be used to 
estimate the limits on the total molecular gas mass, Mn^, 
in HDF850.1. Adopting the relationship: AfH2 = Q^-^co(i-o) 
(q ~ 1 M0(K km s~^ pc^)~^), appropriate for nearby ul- 
traluminous infrared galaxies (Downes & Solomon 1998), 
we can estimate an upper limit to the molecular gas mass 
contained within HDF850.1, under the assumption that its 
redshift is in the range, 3.3 ^ z ^ 5.4. The 3-a limit 
to the CO J=l-0 line luminosity is in the range, L'qq <, 
(3.7 — 8.3) X 10^" K km s~^ pc"^, depending on the assumed 
line width and redshift. By first accounting for a lensing am- 
plification factor of 3 (Dunlop et al. 2004), these limits on 
the CO line luminosity translate directly to a molecular gas 
mass, A'/h2 ^ (1-2 - 2.8) x 10^° Mq. The lensing amplifi- 
cation factor has been calculated for the FIR/submm emis- 
sion region, which we are assuming is co-spatial with the 
low- J CO line emission region, resulting in an equal lensing 
factor. This assumption may not be valid, but without high 
angular resolution observations of both the FIR/submm and 
the low- J CO emission line regions, it remains uncertain. 



3.2 CO and Far-Infrared luminosities 

We can assess whether the CO line luminosity limits 
achieved here, are sufficient to have detected CO line emis- 
sion over the observed redshift interval, by considering the 
CO line luminosity (L'qq) predicted by the estimated far- 
infrared luminosity (Lfir) of HDF850.1 within the context 
of the locally observed Lfir — L'qq relation. In nearby galax- 
ies there exists a well-established correlation between far- 
infrared luminosity and CO line luminosity (e.g. Young & 
Scoville 1991), though it is unclear whether this relation 
truly arises from a direct dependance of star-formation rate 
(as traced by Lfir) on the total molecular gas mass (as 
traced by Lqq)- Furthermore, this relationship appears to 
deviate from a power-law at high far-infrared luminosities 
(Lfir ;a lO''' Lq), which are characteristic of the SMC pop- 
ulation. Despite the uncertainties in this relation, we con- 
verted the estimated far-infrared luminosity in HDF850.1 
to an expected CO J=l-0 line luminosity, in order to de- 
termine if our limits on the CO J=l-0 and CO J=2-l line 
luminosities are sufficiently sensitive for us to have confi- 
dently expected a CO detection. 

Following Neri et al. (2003) and Greve et al. (2005), we 
calculate the far- infrared luminosity for HDF850.1 accord- 
ing to, Lfir ~ 1.9 x 10^^ SgsoimJy] L© (Blain et al. 2002), 
under the assumption of a modified greybody with dust tem- 
perature Td = 40 K, and emissivity index /3 = 1.5, where 
5*850 is the observed 850 flux density. Although various 
measurements of the 850 fim fiux density in HDF850.1 exist 
in the literature, the differences are not significant and we 
adopt the original value presented by Hughes et al. (1998), 
5*850 ~ 7.0±0.5 mjy. This leads to an estimated far-infrared 
luminosity, Lfir ~ (13.3±1.0)m~^ x 10^^ Lq (where m is the 
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Figure 4. a) K-band spectra of HDF850.1 produced from the best 47% of the data, or 12.3 and 13.9 hours of on-source integration time 
included in the left and right polarization spectra, respectively. Channel widths are 8.786 MHz. The black solid-line shows the scaled 
K-band receiver temperature sampled at 30 MHz resolution. Receiver resonance lines appear at various frequencies across the band. 
Grey boxes show the 6 frequency tunings (each ~1.5 GHz wide) adopted to cover the full K-band window. Receiver resonances, such 
as the one at 25.7 GHz, also appear in the HDF850.1 spectra. Another noteworthy feature is the emission line at 19.6 GHz seen in the 
left polarization spectrum but not in the receiver temperature curve. This emission feature is also observed in the calibration spectra 
of 3C295, so we believe it to be another receiver resonance line, b) The 3-cr upper-limits to the CO J=l-0 line luminosity calculated 
from the co-added left and right polarization spectra, assuming CO line widths of 1100, 780 and 460 km s~^. Also shown as the shaded 
region, is the range encompassed by the first and third quartiles of the CO line luminosities in the first 12 SMGs detected in high- J CO 
line emission (see Greve et al. 2005). 



magnification factor due to gravitational lensing, believed to 
be ~3; Dunlop et al. 2004). 

In Figure [S] we compare the estimated limits on the 
ipiR — L'qq parameter space obtained here for HDF85Q.1, 
with the I/FiR — L'qq relation observed in other SMGs and 
various low-redshift galaxy samples. The ultraluminous in- 
frared galaxy (ULIRG) sample observed in CO J=l-0 by 
Solomon et al. (1997) is also included, along with the SLUGS 
sample of Dunne et al. (2Q0Qb) with CO J=l-0 line lumi- 
nosities taken from the literature (Sanders et al. 1985, 1986, 
1991; Young et al. 1995; Casoh et al. 1996; Chini, Kriigel 
& Lemke 1996; Maiohno et al. 1997; Solomon et al. 1997; 
Lavezzi & Dickey 1998). More recent measurements of the 
nuclear CO J=l-0 and CO J=3-2 line emission in a sub- 



set of the SLUGS are presented in Yao et al. (2003). The 
SMG sample consists of those objects in which searches have 
been conducted for high- J CO line emission (Prayer et al. 
1998, 1999; Neri et al. 2003; Genzel et al. 2003; Sheth et 
al. 2004; Greve et al. 2005; Kneib et al. 2005; Tacconi et 
al. 2006), and the one object, SMM J13120-I-4242 at z=3.4, 
for which the CO J=l-0 line has been detected (Hainline et 
al. 2006). The far- infrared luminosities of these objects are 
calculated following the same prescription as that adopted 
for HDF850.1, with submm/mm flux densities taken from 
the literature (Smafl et al. 1997, 1998; Ivison et al. 1998; 
Barger, Cowie, & Sanders 1999; Dey et al. 1999; Cowie, 
Barger, & Kneib 2002; Scott et al. 2002; Chapman et al. 
2003b, 2005; Greve et al. 2004). For the purpose of com- 
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Figure 5. a) The Q-band spectra of HDF850.1 showing only the best of the combined left polarization spectra sampled at 7.813 MHz 
resolution. The right polarization spectra suffered from a severe baseline ripple of unknown origin so are not included here. The redshift 
range (z ~ 3.91 - 4.25) covered by this search for CO J=2-l improves upon the K-band spectra centered at 22.5 GHz, where receiver 
resonance lines made it difficult to place constraints on the presence of CO J=l-0 line emission, b) The luminosity limits on CO J=2-l 
derived from the noise in the left polarization, Q-band spectrum, assuming CO linewidths of 1100, 780 and 460 km . As in Figure|4lb), 



we also show the range defined by the first and third quartiles of the CO line luminosity (median L'^q 
given in Greve et al. (2005). 



X 10^° K km s-l pc2). 



parison with the nearby galaxies detected in CO J=l-0 
, we follow Greve et al. (2005) and assume that for the 
high- J CO lines observed in SMGs, ^co{4-3)/^co(i-o) ~ 
-'^co(3-2)/-^''co(i-o) — -^'co(2-i)/-t'co{i-o) = 1> Correspond- 
ing to optically-thick, thermalized CO emission. These data 
are plotted in Figure |S] with appropriate corrections ap- 
plied to both I/piR, and L'qq to account for magnification by 
gravitational lensing in the 6 SMGs believed to be lensed 
(assuming co-spatial far-infrared and CO emission line re- 
gions). 

Also plotted in Figure [6] is a fit to the Lfir — Leo 
relation derived from a combined sample of ULIRGs and 
SMGs by Greve et al. (2005) where log Leo = (0-62 ± 
0.08) log Lfir + (2.33 ± 0.93). We find a similar relation 
when we fit to the luminosities in the combined low-redshift 
LIRGs/ULIRGs and SMC sample plotted here, log Leo = 
(0.57 ± 0.03) log Lfir + (3.10 ± 0.34). These fits may be 



used to compare the estimated far-infrared luminosity in 
HDF850.1 to the measured CO line luminosity limits. Af- 
ter correcting for amplification by gravitational lensing, 
the estimated far-infrared luminosity in HDF850.1 makes 
it one of the least intrinsically luminous SMGs that has 
been searched for CO line emission, as most in the sam- 
ple are believed to be unlensed (see Greve et al. 2005). 
Adopting the estimate for Lfir in HDF850.1, the fit to the 
Lfir — L'qq relation by Greve et al. (2005) would predict, 

~^ pc^, while the fit presented 



^co 



K km s"^ pc^ 



here would predict, L'qq = 3.4m~ x 10^ 
These values are generally lower than the 3-a CO line lumi- 
nosity limits achieved in our K-band and Q-band spectra. 

We are unable to draw any conclusions as to the validity 
of the photometric redshift technique applied to HDF850.1. 
Although we have removed ~53% of the original K-band 
data due to various spectral baseline irregularities, a further 



8 Wagg et al. 



T ' I 




LpiR [Ld 

Figure 6. The relation between CO line luminosity and far- 
infrared luminosity in various samples of luminous infrared galax- 
ies and AGN. The diamonds are from the CO J=l-0 line 
observations of nearby ULIRGs presented in Solomon et al. 
(1997), while the squares are from the SLUGS sample of nearby 
LIRGs/ULIRGs observed at submm wavelengths by Dunne et al. 
(2000b). The solid circles and short arrows represent the sample 
of high-redshift SMGs which have been observed in high-J CO 
line emission (where both Lpm and L'^q has been corrected for 
lensing when applicable) and the solid triangle is for the one 
SMG detected in CO J=l-0 line emission (Hainline et al. 2006). 
The range in limits to the CO line luminosity in HDF850.1 mea- 
sured here (under the assumption of a redshift, 3.3 ^ 2 ^ 5.4) are 
shown as the long solid arrow, and have been corrected for a lens- 
ing magnification factor of 3. We show the fit presented in Greve 
et al. (2005) to a ULIRG-I-SMG sample {dashed line), along with 
the fit derived here to the low-redshift objects plotted, and only 
those SMGs which are detected in CO line emission (dot-dashed 
line). 

-v/2 decrease in the noise would still not be sufficient to ob- 
tain an L'uQ limit that was significantly inconsistent with 
the expected CO content, given tiie uncertainties and the 
dispersion in tlie estimated gas masses of SMGs with CO 
detections. The CO J=l-0 and CO J— 2-1 line luminosity 
Umits presented here are not of sufficient depth to exclude 
the presence of CO line emission within the redshift interval, 
3.3 ^ z <: 5.4. 



4 CONCLUSIONS 

We present a broadband, GBT spectroscopic search for 
CO J=l-0 and CO J=2-l line emission in the high-redshift 
SMG, HDF850.1 using the K-band (18.0 to 26.5 GHz) and 
Q-band (40.0 to 48.0 GHz) receivers. Although we do not de- 
tect any CO line emission in this object, our constraints on 
the CO line luminosity are approaching that predicted by 
the far-infrared luminosity within the context of the local 
Lfir — L'qq relation. These GBT results are still consistent 
with HDF850.1 lying in the redshift interval, 3.3 ^ 2 ^ 5.4, 
based on our previous rest-frame far-infrared-to-radio pho- 
tometric redshift estimate. 

The GBT has recently been successful in detecting the 
CO J=l-0 line in 3 quasar host galaxies; APM 08279-1-5255 
at z=3.9, PSS J2322-hl944 at z=4.1 and BR 1202-0725 at 



z=4.7 (Riechers et al. 2006), while Hainline et al. (2006) 
present a detection of CO J=l-0 line emission in the SMG 
SMM J13120-f4242 at z=3.4. These 4 objects detected in 
CO J=l-0 line emission with the GBT, were previously 
known to exhibit strong high- J CO line-emission, and there- 
fore to contain large masses of warm molecular gas. With 
these prior high- J CO detections, the redshifts for the molec- 
ular emission-line regions were constrained to ;$100 km s~^ 
enabling a more efficient use (i.e. a narrower frequency 
search) of their available GBT observing time. Given the 
prior uncertainty in both the redshift and CO line inten- 
sity of HDF850.1, our experiment is quite distinct from the 
'tuned' GBT observations described above. This is the first 
broad-bandwidth cm-wavelength search for CO-line emis- 
sion in a high-redshift object (guided by a radio-to-FIR pho- 
tometric redshift), with no previous detections of molecular 
line emission or an optical redshift. 

Considering future possibilities, the gaseous medium 
within SMGs is expected and perhaps already shown to 
be warm and dense, and hence the high-J CO line tran- 
sitions should be more intense than the J=l-0 transition 
(Weiss et al. 2005). At z~2 - 4 (typical of the SMG pop- 
ulation) the high-J (J 2) CO transitions are redshifted 
into the ~70 — 310 GHz atmospheric windows. Thus we 
are optimistic that CO line searches, using broadband mm- 
wavelength receivers on sensitive facilities such as the LMT, 
PdBI or CARMA, will be more successful in obtaining un- 
ambiguous spectroscopic redshifts for the optically obscured 
SMG population of starburst galaxies. 
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